Primordial black holes as a tool for constraining non-Gaussianity 
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Primordial Black Holes (PBH's) can form in the early Universe from the collapse of large density 
fluctuations. Tight observational limits on their abundance constrain the amplitude of the primordial 
fluctuations on very small scales which can not otherwise be constrained, with PBH's only forming 
from the extremely rare large fluctuations. The number of PBH's formed is therefore sensitive to 
small changes in the shape of the tail of the fluctuation distribution, which itself depends on the 
amount of non-Gaussianity present. We study, for the first time, how quadratic and cubic local 
non-Gaussianity of arbitrary size (parameterised by /nl and <?nl respectively) affects the PBH 
abundance and the resulting constraints on the amplitude of the fluctuations on very small scales. 
Intriguingly we find that even non-linearity parameters of order unity have a significant impact on 
the PBH abundance. The sign of the non-Gaussianity is particularly important, with the constraint 
on the allowed fluctuation amplitude tightening by an order of magnitude as /nl changes from just 
—0.5 to 0.5. We find that if PBH's are observed in the future, then regardless of the amplitude of 
the fluctuations, non-negligible negative /nl would be ruled out. Finally we show that <;nl can have 
an even larger effect on the number of PBH's formed than /nl. 



PACS numbers: 95.35.+d 



I. INTRODUCTION 



Primordial Black Holes (PBH's) play a very special 
role in cosmology. They have never been detected but 
this very fact is enough to rule out or at least tightly 
constrain many cosmological paradigms. Convincing the- 
oretical arguments suggest that during radiation domi- 
nation they can form from the collapse of large density 
fluctuations If the density perturbation at horizon 
entry in a given region exceeds a threshold value, of or- 
der unity, then gravity overcomes pressure forces and the 
region collapses to form a PBH with mass of order the 
horizon mass. 

There are tight constraints on the abundance of PBH's 
formed due to their gravitational effects and the con- 
sequences of their evaporation (for recent updates and 
compilations of the constraints see Refs. These 
abundance constraints can be used to constrain the pri- 
mordial power spectrum, and hence models of inflation, 
on scales far smaller than those probed by cosmological 
observations (e.g. Refs. 0-0]). These calculations usu- 
ally assume that the primordial fluctuations are Gaus- 
sian. However since PBH's form from the extremely rare, 
large fluctuations in the tail of the fluctuation distribu- 
tion non-Gaussianity can potentially significantly affect 
the number of PBH's formed. Therefore PBH formation 
probes both the amplitude and the non-Gaussianity of 
the primordial fluctuations on small scales. 

Bullock and Primack Q and Ivanov Q were the first 
to study the effects of non-Gaussianity on PBH for- 
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mation, reaching opposite conclusions on whether non- 
Gaussianity enhances or suppresses the number of PBH's 
formed (see also Ref. [l0[). Refs. [HI [ljj used a non- 
Gaussian probability distribution function (pdf) derived 
from an expansion about the Gaussian pdf [TH, [T3| to 
study PBH formation. However, since PBH's form from 
rare fluctuations in the extreme tails of the probability 
distribution, expansions which are only valid for typical 
size fluctuations can not reliably be used to study PBH 
formation. 



Ref. [HI studied the constraints from PBH formation 
on the primordial curvature perturbation for the special 
cases where it has the form £ = ±(x 2 — (x 2 )), where x 
has a Gaussian distribution, see also Ref. [1 6I | . The minus 
sign is expected from the linear era of the hybrid inflation 
waterfall (see also Ref. (T3] ) , while the positive sign might 
arise if £ is generated after inflation by a curvaton-type 
mechanism. 



In this paper we go beyond this earlier work and cal- 
culate the constraints on the amplitude of the primordial 
curvature fluctuations, £, from black hole formation for 
both the quadratic and cubic local non-Gaussianity mod- 
els (parameterised by /nl and c/nl respectively). In the 
process we calculate the probability distribution function 
of the curvature perturbation for these models. Our re- 
sults are valid for arbitrary values of these non-linearity 
parameters, and we recover the known limiting results 
for very small or large non-Gaussianity. In Sec. [II] we re- 
view the calculation of the PBH abundance constraints in 
the standard case of Gaussian fluctuations, before calcu- 
lating the constraints for quadratic and cubic local non- 
Gaussianity in Sec. Mil and IIVI respectively. We conclude 
with discussion in Sec. fV] 
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II. PRIMORDIAL BLACK HOLE FORMATION 
CONSTRAINTS 

The condition for collapse to form a PBH is tradition- 
ally stated in terms of the smoothed density contrast at 
horizon crossing, <5i lor (i?). A fluctuation on a scale R 
will collapse to form a PBH, with mass Mpbh roughly 
equal to the horizon mass, if <5hor(-R) > 5 C ~ 0(1) [l[ 1 . 
If the initial perturbations have a Gaussian distribution 
then the probability distribution of the smoothed density 
contrast is given by (e.g. Ref. [lH): 



P(S h or(R)) = 



1 



■ exp 



where cr(R) is the mass variance 



, (1) 



a\R) = J W z (kR)T s (k,t) T , (2) 

while Vs(k, t) is the power spectrum of the (unsmoothed) 
density contrast 



.2/n\ _ / TT/2/- — — " - aK 



p s (k,t) = ^(\S k \ 2 ), 



(3) 



and W(kR) is the Fourier transform of the window func- 
tion used to smooth the density contrast. 
The initial PBH mass fraction 



PBH ) 



_ Ppbh(^pbh) 

Ptot 



(4) 



is equal to the fraction of the energy density of the Uni- 
verse contained in regions dense enough to form PBH's 
which is given, as in Press-Schechter theory (20| . by 



a wide range of constraints on the PBH abundance, from 
their various gravitational effects and the consequences of 
their evaporation, which apply over different mass ranges. 
These constraints are mass dependent and lie in the range 
/3(M PBH ) < lfr 20 -lCr 5 @,[|. The power of these PBH 
abundance constraints is apparent when we consider the 
resulting constraints on cxh or (i?) which are in the range 
Phor(R) / $c < 0.1 — 0.2. In other words a small change in 
0hor{R)/8c leads to a huge change in j3. Finally to im- 
pose constraints on the power spectrum of the curvature 
perturbation the transfer function which relates 5(k, t) to 
the primordial curvature perturbation is calculated (e.g. 
Refs. Hip). 

For an interesting number of PBH's to form the power 
spectrum on small scales must be several orders of mag- 
nitude larger than on a cosmological scales. This is pos- 
sible in models, such as the running mass model |22h27| . 
where the power spectrum increases monotonically with 
increasing wave- number [f| Q. Another possibility is a 
peak in the primordial power spectrum, due to a phase 
transition during inflation or features in the inflation po- 
tential (see Ref. [28| and references therein). 

To assess the affects of non-Gaussianity on the PBH 
constraints on the curvature perturbations it is sufficient 
to follow Ref. [lH and work directly with the curvature 
perturbation, with the PBH abundance being given by 



^(C)dC- 



For Gaussian fluctuations 
1 



P(0 



2na 



cxp 



2a 2 J 



(7) 



(8) 



and hence 



/3(M 



PBH J 



<5c 



P(5 hor (R))d5 hm (R) 



(5) 



where we have followed the usual Press-Schechter prac- 
tice of multiplying by a factor of 2, so that all the mass 
in the Universe is accounted for. The PBH initial mass 
fraction is then related to the mass variance by 



/3(M] 



PBH J 



erfc 



27TCThor(-R) J 5c 

5c 



exp 



ghorCg) 



dS hor (R) . 



V2a hm (R) 



[1 = 2 erfc 



V2cr 



(9) 



where Cc is the threshold for PBH formation. The vari- 
ance of the probability distribution is related to the power 
spectrum of the curvature perturbation, 



2tt 2 



(10) 



(6) 



The constraints on the PBH initial mass fraction, 
/3(Mpbh), can therefore be translated into constraints on 
the mass variance by inverting this expression. There are 



1 It was previously thought that there was an upper limit on the 
size of fluctuations which form PBH's, with larger fluctuations 
forming a separate closed universe. Kopp et al. have recently 
shown that this is in fact not the case. 



by a 2 « "Pf . Here and subsequently for compactness we 
drop the explicit scale dependence of ft and a and the 
subscript 'hor' indicating that a is to be evaluated at 
horizon crossing of the mass variance, i.e. a = ^^(R). 

The constraints on a obtained via this method differ 
from those obtained from the full calculation involving 
the smoothed density contrast by 0(10%): for /3 = 10~ 20 
the full calculation gives V^ 2 = 0.12 0], while using 
Eq. © gives a = 0.11. We take £ c = 1 for defmiteness, 
however any variation in the threshold for collapse affects 
the constraints on a in the Gaussian and non-Gaussian 
cases in the same way (since it is in fact the combination 
u/C,c that is constrained). We consider the upper and 
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lower values of the constraints on /3, 10~ 5 and 10 -20 , 
and henceforth drop the explicit dependence on the PBH 
mass. 

Because PBH's form from the rare large fluctua- 
tions in the tail of the distribution, this corresponds to 
Q c /{\/2g) 3> 1 which allows a useful analytic approxi- 
mation for a(f3) to be found. Using the large x limit of 
crfc(x), Eq. ([S]) can be rewritten as 



exp 



, 2 
2a 2 



(11) 



The non-Gaussian probability distribution function 
(pdf) for C, Png(C), can b e found by making a formal 
change of variables 



p NG (C)dC = ^ 



PaiK 1 )^, (15) 



where Pg has a Gaussian distribution and the sum is 
over the n solutions, of the equation h((a) — C 

For quartic local non-Gaussianity solving Eq. (fTB"]) for Cg 
gives two solutions (which we denote with '±'): 



and hence 



Cc 



2 Ml/0) 



(12) 



Up to now we have been concentrating on the large tail 
limit of a Gaussian distribution. We will now see that 
the consequence of allowing even a small amount of non- 
Gaussianity in the distribution can be dramatic in terms 
of the constraints from PBH formation - enough to po- 
tentially rule out certain values of the non-linearity pa- 
rameters. Wc begin in Scc. lIIII bv considering the impact 
of adding quadratic local non-Gaussianity. 



III. QUADRATIC LOCAL NON-GAUSSIANITY 

(/nl) 

We take the model of local non-Gaussianity to be 



C = Co + g/NL(CG - O = MCg), 



(13) 



where Cg is Gaussian with variance a, i.e. (C G ) = a 2 . 
The constant term is subtracted from C such that (C) =0, 
as required by the definition of the curvature perturba- 
tion. Note that we are not assuming single-field infla- 
tion, but the less restrictive assumption that only a single 
field direction generates C, known as single-source infla- 
tion [29j . and furthermore this degree of freedom might 
be different from the one which generates the primor- 
dial temperature perturbation on the much larger CMB 
scales, this is for example the case in [l5j]. The variance 
of C is given by [3(| 



(C 2 ) 



3./': 



NL 



er 4 ln(fcL), 



(14) 



where the cut-off scale L ~ 1/H is of order the hori- 
zon scale, the second term is a one-loop correction which 
dominates in the large /nl limit, k is the scale of interest 
and ln(fc-L) is typically of order unity (see e.g. Ref. f3lj]). 
We have implicitly assumed here, for simplicity, that the 
primordial power spectrum is close to scale-invariant on 
the scales of interest, as is the case for models where the 
power spectrum varies monotonically with wavenumber. 



6/nl 
and hence [lj, |3 

^ng(C)c1C = 



(,16) 



3/nlct 2 _|_ £ 



where 



e± = exp 



1 (h±\Q 



(e-+e_) , 
(17) 

(18) 



The log of the non-Gaussian pdf is shown in the upper 
panels of Fig. [1] for a = 0.1 and / NL = 0, ±2, ±3.5 and 
±5. For positive /nl, as /nl is increased the amplitude 
of the large £ tail of the pdf increases in amplitude. There 
is a minimum value of C, Cmin = Ciim, 







lim 



12/ 



NL 



1 



36/nl^ 
25 



(19) 



at which the pdf diverges, since d/i _1 (C)/dC tends to in- 
finity as C — > Ciim- Initially the peak in the pdf increases 
in amplitude and moves to negative £. As /nl is in- 
creased further the pdf becomes monotonic, increasing 
continuously with decreasing £ down to Cmin- For nega- 
tive /nl the pdf is the mirror image of that for positive 
/nl, and in this case there is a maximum possible value 

of C, Cmax Clim- 

The initial PBH mass fraction is given by 



P = !++!-, 



where 



27R7J1 



_|_ 12/nl ^ 3/nlo- 2 _|_ £j 



(20) 



=e±dC-(21) 



The lower limit on the integral is Cc — 1, the threshold 
value above which a black hole is formed. The upper 
limit, Cmax, is the maximum possible value of £. For 
/nl > 0, Cmax = oo, while for /nl < 0, Cmax = Ciim as 
discussed above, with Ciim given by Eq. (|T9"1) . 
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FIG. 1: The log of the non-Gaussian probability distribution, 
top row /nl and bottom row jnl- We have fixed a — 0.1 
throughout and the solid line shows the Gaussian pdf. In the 
top left panel the dotted, short-dashed and long-dashed lines 
correspond to /nl = 2, 3.5 and 5 respectively, while in the top 
right panel they correspond to /nl = —2, —3.5 and —5. In the 
bottom left panel the dotted, short-dashed and long-dashed 
lines correspond to jnl = 10, 20 and 30 respectively, while in 
the bottom right panel they correspond to <?nl = —10, —20 
and —30. 
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FIG. 2: The constraints on the square root of the power spec- 

1/2 

trum of the curvature perturbation, , for initial PBH 
abundances /3 = 10~ 5 and 10 -20 (lower and upper lines re- 
spectively) for the quadratic local non-Gaussianity model as 
a function of /nl. 



The initial PBH mass fraction is most easily calculated 
by making a transformation to a new variable y, 



h ± H0 



FIG. 3: The constraints on the square root of the power spec- 
trum of the curvature perturbation, , as in Fig. [21 for 
the quadratic local non-Gaussianity model as a function of 
the fraction of the power spectrum which is non-Gaussian, 
4sgn(/)/ 2 a 2 /(l + 4/ 2 ( 7 2 ), where / = 3/ NL /5. 



which has unit variance so that 
./nl > 0: 







/nl < 0: 
1 




dy + f~ e-*'/ 2 dy j (23) 

J — DO J 

(24) 



e~» ' /2 dy, 



(25) 



where are the values of y corresponding to the thresh- 
old for PBH formation, £ c : 



y± 



(26) 



(22) 



For /nl >0, y% > 0, y c _ < 0, and = y% + y c _ = 

—3/ (5/nl)- Consequently the first integral in the expres- 
sion for j3, Eq. (f23|) . which corresponds to the positive 
branch, gives the dominant contribution to j3. However 
in the limit of very large /nl j | y+ \ ~ \ V- \ tends to zero and 
the positive and negative branches contribute equally to 
/?. 

The constraints on the square root of the power spec- 

1/2 

trum of the curvature perturbation, , which arise 
from the tightest and weakest constraints on the initial 
PBH mass fraction, /3 < 10~ 20 and 10 -5 respectively, are 
shown in Figs. [2] and [3] In Fig. [2] we plot the constraints 
as a function of /nl for small /nl, while in Fig. [3] we 
plot the constraints as a function of the fraction of the 
power spectrum which is non-Gaussian i.e the ratio of the 
second non-Gaussian term in the expression for (£ 2 ) in 
Eq. (Q3J) to the full expression. The limit that this ratio 
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is unity corresponds to a purely non-Gaussian C- We now 
discuss how the change in the PBH constraints depends 
on the amount of non-Gaussianity. 



A. Very small |/nl| 

Expanding the pdf for £, Eq. (|T7| . to second order in 
/nl wc find 



Png(C)=^g(C) 



(3/ NL C) 2 fC 4 „C 



50 



- 11 





-') 


3/nlC 








23- 


4)} 



, (27) 



which agrees with the non-Gaussian pdf in Ref. [13[ to 
linear order in /nl- The second order term has an extra 
l/er 2 term, and hence the above expansion can only be 
expected to be accurate for /nl <C <t 2 /C $ - Since Cc ~ 1 
and a ~ 0.1 in the Gaussian case, this means that this 
expansion is only valid when applied to PBH formation 
if /nl < 0.01. 

As can be seen in Fig. [2] even a small level of non- 
Gaussianity has a significant impact on the constraints 
on a. This is because PBH's form from the fluctuations 
in the extreme tail of the distribution (e.g. for a ~ 0.2 
and (^c = 1 they are a 5tr fluctuation) and it is in this 
regime that even a small skewness is important. The 
strong asymmetry between positive and negative /nl is 
because for /nl > overdensities are enhanced, while 
for /nl < the overdensity from a positive C is partially 
canceled by the Gaussian squared term, thereby a has to 
become large in order for PBH formation to be possible 
at all. 



B. Intermediate |/nl| 

In the regime where < /nl"" 1 (which in practice 
corresponds to < /nl < 10 since a ~ 0.1), we have 
CT 2 ^ However /nlCc ^ 1 and the expression for 

/i _1 (C), Eq. (|TB)) . simplifies substantially leading to 



a 

Cc 



3/nlCc 



1/2 



21n(l//3) 



(28) 



The square root term is the result in the Gaussian case, 
Eq. (|T2l) . Hence the constraint on V( is tightened by a 
factor of /nl compared to the Gaussian result. 

As can be clearly seen from Fig. [5J the constraints are 
very asymmetric under a change of sign of /nl, with 
the constraints becoming very rapidly much weaker for 
/nl < 0, wc discuss the case of negative /nl in the next 
section. 



C. Large |/ N l| 

In the case of a pure, positive non-Gaussianity the con- 
straints on become a lot tighter. Since there is a de- 
generacy between /nl and a in this case, £ can be taken 
to be given by 



C = ± (Co " ° 2 ) = HCg) < a 2 



(29) 



(c.f. Ref. (3) and hence (C 2 ) ^ 4cr 4 . 

We first study the + case. Performing a similar calcu- 
lation to the more general case with a linear term, the 
PBH initial mass fraction is given by 



and hence 111 



-1/72 



dy, 



1/2 



2(7 2 



Cc 



ln(l//3)' 



(30) 



(31) 



i.e. the constraint on a is approximately the square of 
the constraint in the Gaussian case, and is hence a lot 
tighter. 

The case of a pure, negative chi-squared distribution is 
very different from the positive case. Using the transform 
h^ 1 this leads to 



Cc 



2^v^^ 



dC: 



(32) 



which using a further change of variables we transform 
to 







1 



-y/2 



dy, 



<T 2 -Cc 



(33) 



which leads to the relationship between a 2 

a 2 = Cc (1 + 8™ 2 /? 2 Cc + 0(/? 4 )) * C< 



Vl 12 and (3 



(34) 



The constraint on a is very weakly dependent on the 
limit on /3, confirming the behaviour seen in Fig. [2] 
where the constraints for j3 = 10~ 5 and 10 -20 merge for 
/nl £ -0.5. For / NL < 0, once \ f NL \a > 0.4 the pdf for 
C increases monotonically with increasing C, before di- 
verging at Cmax = Ciim- If Cmax < Cc then the number of 
PBH's formed is identically zero, while if Cmax > Cc it is 
extremely sensitive to the precise value of a. Therefore in 
this regime, unless there is extreme fine-tuning of a, the 
number of PBH's formed will cither be completely negli- 
gible or so large that the PBH abundance constraints are 
violated by many orders of magnitude. Although one can 
formally produce any given value of (3 with sufficient fine 
tuning of a, in a realistic model the non-Gaussianity will 
lead to small spatial variations of a in different patches 
(e.g. due to a small cubic term ) [33[, which would prob- 
ably rule this model out if one performed a more detailed 
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calculation. Hence we conclude that a future detection 
of PBH's would effectively rule out a negative /nl unless 
it has a tiny value, i.e. from Fig. [21 /nl ^ —0.5 would be 
ruled out regardless of the value of a. 

Having looked at the case of adding a quadratic type of 
local non-Gaussianity, we now consider the case of adding 
a cubic type to see what new constraints this may impose. 



IV. CUBIC NON-GAUSSIANITY (# N l) 

The model of local non-Gaussianity with a cubic term 
(but assuming that /nl = 0) is defined by 



C = Cg+.9Cg = MCg 



25 5NL ' 



(35) 



where we have introduced the definition of g in order to 
reduce the numerical factors which will appear in many 
expressions in this section. The variance of C for this 
model is given by 

(C 2 > = o- 2 (1 + 6g<r 2 hi(kL) + 27. 9 V ln(fcL) 2 ) , (36) 

where the second term is a one-loop contribution and the 
second term a two- loop contribution (30| , which nonethe- 
less dominates in the limit of large g, and ln(fcL) is again 
of order unity. 

For (?nl > the cubic equation C = Cg + ffCc = MCg) 
has one real solution for all £: 



2 l/3 



g z (C + \K 2 + 



27g 



-1/3 



+ 



1 



9 Z ic+\ e + 



273 J 



(37) 



The PBH initial mass fraction is then given by 
g > 0: 



13 = 



1 



,-y 2 /? 



dy 



erfc [ % 
,V2 



(38) 



where y c = h~ 1 (( c )/a, with /i _1 (C c ) given by Eq. (j3"7) . 

In the case of negative <7nl, the cubic function has 
a local maximum for positive Cg, with a peak value at 
C = Ct where 



Ct 



3V3 



9 



(39) 



Hence the cubic polynomial has three real roots if |C| < 
Ct, and otherwise only one real root. The transition be- 
tween one and three real roots for C = Cc occurs when 
g = g t where 



27C 2 



-0.15. 



(40) 



The expression for the PBH initial mass fraction there- 
fore has different forms depending on whether g is greater 
or less than g t . 



If Ct < Cc, cquivalently g < gt, then 
9 < 9t ■ 



1 



/2tt 

where y c = /i _1 (Cc)/c and 
2 i/3 



cxp (-y 2 /2) dy = erfc ( ^ 



3(- 5 ) 2 /3 
1 

"2l/3(- fl )l/3 



Cc + vc 2 -^ 



^/2 



-1/3 



(41) 



Cc + \/C 2 -Ct 



1/3 



(42) 



If Ct > Cc, cquivalently g t < g < 0, there are three 
roots ftf^Cc) < < h^iCc) < h^(Cc), given by 



K\Q = 

K\Q = 

K\Q = 

where 



Vsi-g) 1 / 2 
1 

\/3(-5) 1/2 

1 

V3(-s) 1/2 



cos {9/3), (43) 
cos (6>/3) - V3sin(6»/3)j (44) 
cos(6»/3) + V3sin(6»/3)j (45) 



= at an 



(C?-£) 



2)1/2 



Cc 



(46) 



It follows that 
gt < g<0: 

P = 



1 ( f Vl 

= (j cxp(-y 2 /2)dj/ 
f 3 exp {-y 2 /2) dy ) , 



(47) 



erfc 



V2 



cxp(V/2)dy, (48) 



U-2 



where j/j = h i 1 {C tC )/a, with h i 1 {C,c) given by Eqs. (|4"3jl- 

We calculate the non-Gaussian pdf using the procedure 
described in Sec. IIII1 The log of the pdf is shown in the 
lower panels of Fig. Q] for a = 0.1 and <?nl = 0, ±10, ±20 
and ±30. For positive <7nl, as <7nl is increased the large 
C tail of the pdf increases in amplitude, however (un- 
like the /nl case) the body of the pdf does not deviate 
significantly from the Gaussian pdf. This suggests that 
PBH formation is potentially a more sensitive probe of 
positive cubic local non-Gaussianity than structure for- 
mation. Negative c/nl is very different from positive <?nl- 
In particular there is a divergence at C = Ct which arises 
from the d/i -1 (C)/dC factor in the pdf. From Eq. (|3"5)) we 
see that 



dC 



1 + 3gC 2 G 



(49) 
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Snl 

FIG. 4: The constraints on the square root of the power spec- 
trum of the curvature perturbation, P^ 2 , as in Fig.O for the 
cubic local non-Gaussianity model for — 10 < <?nl < 10. 
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shown in Figs. [4] and [5l In Fig. [4] we plot the constraints 
as a function of (?nl for small <?nl, while in Fig. [5] we plot 
the constraints as a function of the fraction of the power 
spectrum which is non-Gaussian, i.e. the ratio of sum 
of the second and third non-Gaussian terms in Eq. (|36[) 
for (C 2 ) to the full expression. We now discuss how the 
change in the PBH constraints depend on the amount of 
non-Gaussianity. 

A. Small |<7nl| 

The constraints on the power spectrum are highly 
asymmetric between positive and negative (?nl- This is 
because for <?nl > an overdensity in the linear £ regime 
will be boosted by the cubic term, especially strongly in 
the tail of the distribution and hence the constraint is 
tightened. However for mildly negative (?nl, the oppo- 
site is the case; the two terms tend to cancel each other 
out and hence the constraints on the power spectrum 
weaken dramatically in this regime. For very small neg- 
ative <7nl the 2nd term in the expression for /3, Eq. (|47p . 
the integral of the Gaussian distribution from yi to 2/3 
dominates. As g — s- <? t = — 0.15 from above 2/3 — ])2 — >• 
so that this term decreases rapidly and the constraint on 
the power spectrum rapidly becomes weaker. Only when 
<7t — 9 "C 10~ 80 does this term become smaller than the 
first erfc(|j/i|/\/2) term though, and then the value of /3 
matches smoothly onto that in the g < g t regime where 
there is only one real root. As g is decreased below g t 
the one real root, ft. _1 (Cc) given by Eq. (|42|) . becomes 
less negative and the constraint on the power spectrum 
rapidly tightens again. 

B. Intermediate |<?nl| 



FIG. 5: The constraints on the square root of the power 
spectrum of the curvature perturbation, V^ 2 , as in Fig. [3] 
for the cubic local non-Gaussianity model as a function of 
the fraction of the power spectrum that is non-Gaussian, 
sga(g)(6ga 2 + 27 5 2 a 4 )/(l + 6ga 2 + 27g 2 a 4 ). 



This diverges when 1/Cg = 3(— 9) which from Eq. ([55)1 
corresponds to £ = £t ■ Expanding around this point with 
£ = £{ — 6, S <C 1, we find after a little algebra that to 
leading order 



dC 



3-1/4 



2(-5)V4(SVa 



(50) 



This divergence is fairly weak however and the PBH ini- 
tial mass fraction (3 does not diverge, a result that is 
confirmed analytically as well as numerically. 

The constraints on the square root of the power spec- 

1/2 

trum of the curvature perturbation, , which arise 
from the tightest and weakest constraints on the initial 
PBH mass fraction, /3 < 10 -20 and 10 -5 respectively, are 



For (?nl > in the intermediate regime, where the non- 
Gaussianity is small enough that 7^ ~ er 2 , i.e. go 1 ^ 1, 
but large enough to satisfy gQ S> 1, (hence valid for 



1 < g < <t~ 2 ), the expression for /i _1 (Cc) in Eq. (|37j) 
simplifies significantly to /i _1 (C c ) — (Cc/s) 2 ^ 3 - Using the 
leading asymptotic expansion for the Gaussian pdf this 
leads to 



a 



1 

5Cf 



1/3 



1 



21n(l//?)- 



(51) 



The term in the square root is the result in the Gaussian 
case, hence the constraint on a is tightened by a factor 
of .g -1 / 3 . For comparison, the result for intermediate 
quadratic non-Gaussianity is given by Eq. (|28[) . 

To compare the relative size of these changes relative to 
the Gaussian case, consider a 10% non-Gaussian correc- 
tion to £, and for concreteness assume that a = 10~ 2 . 
Then go 1 = 0.1 implies g = 10 3 and the constraint 
on a tightens by a factor of 10. For a quadratic non- 
Gaussianity, 3/nlc/5 = 0.1 implies 3/nl/5 = 10 and 
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hence the constraint only tightens by approximately a 
factor of 3. However, if instead of considering a fixed ra- 
tio of non-Gaussian to Gaussian terms, one considers the 
non-linearity parameters to have a fixed amplitude then 
the constraints on a are tightening by the cube root of 
<?, but only by the square root of / = 3/nl/5. 



C. Large \gai\ 

In the limit of very large |<7nl|, C ^Cg ano - the con ~ 
straints don't depend on the sign of <7nl ■ This is because 
the Gaussian pdf is invariant under a change of sign in 
Qg and this is equivalent to changing the sign of <?nl, 
but only in the case that the linear term is absent. As 
can be seen in Fig. [5] the symmetry between positive and 
negative <?nl only occurs once the modulus of the non- 
Gaussian fraction of the power spectrum becomes close 
to one, which corresponds to very large values of |<?nl|- 
In the limit of very large |<7nl| the constraints on be- 
come significantly tighter than in the case of a large and 
positive /nl- The bound in this case becomes approxi- 
mately 



,1/2 



5a J 



5Cc 



(21n(l//3))3 



(52) 



which is approximately the cube of the bound in the 
Gaussian case and hence much more stringent (and also 
tighter compared to the case of a large and positive /nl, 
as discussed in Sec. MI Cl) . 



V. DISCUSSION 

PBH formation probes the extreme tail of the proba- 
bility distribution function of the primordial fluctuations. 
This is the region of the pdf which is most sensitive to 
the effects of any non-Gaussianity that may be present. 
We have, for the first time, calculated joint constraints 
on the amplitude and non-Gaussianity of the primordial 
perturbations, for arbitrarily large local non-Gaussianity. 
We have studied both quadratic and cubic local non- 
Gaussianity, parameterised by /nl and <?nl respectively 
On the scales associated with the cosmic microwave back- 
ground and large scale structure, the constraints on pri- 
mordial non-Gaussianity are approximately |/nl| ^5 10 2 
34 1 and |</nl| IS 10 6 [35| . In contrast we have shown 
that on much smaller scales non-linearity parameters of 
order unity can have a significant effect on the number 
of PBH's formed. This is because the non-linearity pa- 
rameters have a larger effect on the tails of the fluctu- 
ation distribution than on the more moderate fluctua- 
tions probed by cosmological observations. We expect 
most other forms of non-Gaussianity to also have a sig- 
nificant effect on PBH production, since in general non- 
Gaussianity generates a skewness which affects the tails 
of the pdf. 



The signs of the non-linearity parameters are partic- 
ularly important. If positive they always make the con- 
straints tighter by acting in the same direction as the lin- 
ear contribution to £. A negative quadratic term tends to 
cancel the effect of the linear term, thereby reducing the 
abundance of large PBH forming fluctuations. The con- 
straints on the amplitude of the power spectrum therefore 
become much weaker, of order unity for /nl ^ — 0.5. In 
practice this means that the amplitude of fluctuations 
will either be too small to form any PBH's at all, or so 
large that almost every horizon region collapses to form 
a PBH, which is already observationally ruled out. We 
hence conclude that a future detection of PBH's would 
rule out a negative value of /nl unless its value is tiny, 
|/nl| <C 1- The case of negative <?nl is different. For 
.9nl — — 1 the constraints are weakened as in the neg- 
ative /nl case. However as <?nl becomes more negative 
the constraints quickly become tighter again. In the limit 
of very large <?nl the constraints are independent of its 
sign and very tight, approximately the cube of the con- 
straints in the Gaussian case. We have also studied and 
plotted the probability distribution functions, showing 
that although the pdfs can diverge, all physical quan- 
tities, such as the PBH abundance, remain finite. The 
PBH constraints have previously been calculated for a 
pure x 2 pdf [HI, EH and for quadratic non-Gaussianity 
in the limits that |/nl| *C 1 and the linear term dom- 
inates [ll|, EH- We have shown that we recover these 
limiting cases, however, of particular significance is the 
fact that our calculations are valid for arbitrarily large 
quadratic and cubic local non-Gaussianity. 

A bispectrum in the squeezed limit is in general gen- 
erated by all single field models of inflation, with an 
amplitude which is related to the spectral index by 
/nl = -5(n s - 1)/12 [H, H3- Although this value is 
too small to be seen on CMB scales, it might be impor- 
tant on PBH scales, since firstly the spectral index might 
be larger as the slow-roll parameters potentially become 
larger towards the end of inflation and secondly since the 
constraints are sensitive to smaller values of /nl- 

An important issue, which goes beyond the scope 
of this work, is the calculation of the secondary non- 
Gaussianities generated through the effects of gravity be- 
ing non-linear and through horizon re-entry after infla- 
tion, during which time C is no longer conserved. These 
calculations have been carried out on CMB scales and 
these effects generally cause an order of unity change to 
the non-linearit y p arameters (although the effect is scale 
dependent) [38l - l42l | . Such small values of the non-linear 
parameters can have a significant effect on the number of 
PBH's formed, therefore it would be interesting to carry 
out a similar analysis valid for the much smaller scales 
on which PBH's form. If these effects generically lead to 
/nl ~ — 1 then this would suggest that PBH's are un- 
likely to have formed, unless inflation generated a larger 
and positive primordial /nl on the same scales. The 
current calculation could also be extended by allowing 
for simultaneous non-zero values of /nl and <?nl or by 
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studying the effects of higher order non-linearity param- 
eters. 

The possibility of PBH formation constrains both the 
amplitude and degree of non-Gaussianity associated with 
the primordial density perturbations over a wide range 
of scales, smaller than those probed by cosmological ob- 
servations. In this paper we have shown that even rel- 
atively small values of the non-linearity parameters, of 
order unity, can have a significant effect on the PBH 
constraints on the amplitude of the primordial pertur- 
bations. Therefore non-Gaussianity should be taken into 
account when calculating PBH constraints on inflation 
models. We have also shown that the observation of 



PBH's would rule out non- negligible negative /nl, re- 
cmphasizing the constraining power of PBH's. 
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